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Abstract 

Form factors parameterizing the semileptonic decay Be — > Dsl^l~ {I = /i, r) are calculated using the frame work of Ward 
Identities. These form factors are then used to calculate the physical observables like branching ratio and helicity fractions of 
final state D* meson in these decay modes. The analysis is then extended to the the universal extra dimension (UED) model 
where the dependence of above mentioned physical variables to the compactification radius R, the only unknown parameter 
in UED model, is studied. It is shown that the helicity fractions of D* are quite sensitive to the UED model especially when 
have muons as the final state lepton. Therefore, these can serve as a useful tool to establish new physics predicted by the UED 
model. 
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I. INTRODUCTION 



Living in the LHC era, it is hoped to either verify the Standard Model (SM) or to explore the properties of more 
accurate underlying theory that describes the theory of weak scale. Flavor Changing Neutral Current (FCNC) decays 
of i?-meson are an important tool to investigate the structure of weak interactions and also provide us a frame work 
to look for the physics beyond the Standard Model (SM). This lies in the fact that FCNC decays are not allowed at 
tree level in the SM and occur only at the loop level [l|-|3| and makes them quite sensitive to possible small corrections 
that may be result of any modification to the SM, or from the new interactions. This gives us solid reason to study 
these decays both theoretically and experimentally. 

Since the CLEO observations of the rare radiative 6 — > 57 transition 4], there have been intensive studies on rare 
semileptonic, radiative and leptonic decays of Bu.d,s mesons induced by FCNC transitions of 6 — >■ s, d Q. The study 
will be even more complete if one consider the similar decays of the charmed B mesons (Be). 

The charmed Be meson is a bound state of two heavy quarks, bottom b and charm c, and was first observed in 
1998 at Tevatron in Fermilab Because of two heavy quarks, the Be mesons are rich in phenomenology compared 
to the other B mesons. At the Large Hadron Collider (LHC) the expected number of events for the production 
of Be meson are about 10* — 10^° per year 0, @] which is a reasonable number to work on the phenomenology 
of the Be meson. In literature, some of the possible radiative and semileptonic exclusive decays of Be mesons like 
Be {p,K*,Dl,Bl)j,Be Ivj ,Be ^ B*J+l-,Be ^ 0^1' Be ^ Dl^l+l' and Be Dl J+l- have been studied 
using the frame work of relativistic constituent quark model 0|, QCD Sum Rules and the Light Cone Sum Rules 
The focus of the present work is the study of exclusive Be D*l^l^ decay. 

While working on the exclusive S-meson decays the main job is to calculate the form factors which are the non 
perturbative quantities and are the scalar functions of the square of momentum transfer. In literature the form factors 
for Be D*l^l~ decay were calculated using different approaches, such as light front constituent quark models and 
a relativistic quark model 0, [ll|. In this work we calculate the form factors for the above mentioned decay in a 
model independent way through Ward identities, which was earlier applied to i? — > p, 7 12, 13] and B — > Ki decays 



14j . This approach enables us to make a clear separation between the pole and non pole type contributions, the 
former is known in terms of a universal function ^_l(<;^) = g-\-{q^) which is introduced in the Large Energy Effective 
Theory (LEFT) of heavy to light transition form factors [15]. The residue of the pole is then determined in a self 
consistent way in terms of g+(0) which will give information about the couplings of i?*(l^) and i?*^(l+) with BeD* 
channel. The above mentioned coupling arises at lower pole masses because the higher pole masses of Be meson do 
not contribute for the decay Be — D*l~^l~ . The form factors are then determine in terms of a known parameter g+{0) 
and the pole masses of the particles involved, which will then be used to calculate different physical observables like 
the branching ratio and the helicity fractions of D* for these decays. 

At the quark level the semileptonic decay Be — ^ D'*l^l^ is governed by the FCNC transition b — sl~^l~ , therefore 
it is an important candidate to look for physics in and beyond the SM. Many investigations for the physics beyond 
the SM are now being performed in various areas of particle physics which are expected to get the direct or indirect 
evidence at high energy colliders such as LHC. During the last couple of years there have been an increased interest in 
models with extra dimensions, since they solve the hierarchy problem and they can provide the unified framework of 
gravity and other interactions together with a connection to the string theory [l6j . Among them the special role plays 
the one with universal extra dimensions (UED) as in this model all SM fields are allowed to propagate in available all 
dimensions. The economy of UED model is that there is only one additional parameter to that of SM which is the 
radius R of the compactified extra dimension. Now above the compactification scale 1 /R a given UED model becomes 
a higher dimensional field theory whose equivalent description in four dimensions consists of SM fields and the towers 
of KK modes having no partner in the SM. A simplest model of this type was proposed by Appelquist, Cheng and 
Dobrescu (ACD) 17[. In this model all the masses of the KK particles and their interactions with SM particles and 
also among themselves are described in terms of the inverse of compactification radius R and the parameters of the 



SM Q. 



The most important property of ACD model is the conservation of parity which implies the absence of tree level 
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contributions of KK states to the low energy processes taking place at scale /i << 1/R. This brings interest towards 
the FCNC transitions, 6 — > s as mentioned earlier that these transitions occur at loop level in SM and hence the one 
loop contribution due to KK modes to them could in principly be important. These processes are used to constrain 
the mass and couplings of the KK states, i.e, the compactification radius \/R 3 3- 

Buras et al. have computed the effective Hamiltonian of several FCNC processes in ACD model, particularly in h 
sector, namely Bg d mixing and & — >■ s transition such as & — >■ 57 and b — ^ sl~^l~ decay |18l |. The implications of physics 
with UED are examined with data from Tevatron experiments and the bounds on the inverse of compactification radius 
arc found to be 1/R > 250 — 300 GeV [2^. There exists some studies in the literature on different B to light meson 
decays in ACD model, where the dependence of different physical observables like branching ratio, forward-backward 
asymmetry, lepton polarization asymmetry and the helicity fractions of final state mesons on is examined 20l-l22 1 . 

In this work we will study the branching ratio and helicity fractions of D* meson in Be Dtl~^l~ decay both in 
the SM and ACD model using the framework of i? — > {K* , Ki)l^l^ decays described in refs. [2l|,|22|. The paper is 
organized as follows. In Sec. II we present the effective Hamiltonian for the decay Be D'*l^l^ . Section III contains 
the definitions as well as the detailed calculation of the form factors using Ward Identities. In Sec. IV we present 
the basic formulas for physical observables like decay rate and helicity fractions of D* meson where as the numerical 
analysis of these observables will be given in Section V. Section VI gives the summary of the results. 



II. EFFECTIVE HAMILTONIAN AND MATRIX ELEMENTS 

At quark level, the semileptonic decay Be — 5> D*l~^l~ is governed by the transition b — >■ sl~^l~ for which the effective 
Hamiltonian can be written as 



H, 



eff 



AGf 



VtbV* 
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C.(/i)0, 



(1) 



where Oi(/i) (z = 1,...,10) are the four quark operators and Ci{^) are the corresponding Wilson coefhcients at the 
energy scale n [2^ which was usually take to be the &-quark mass (rrib). The theoretical uncertainties related to the 
renormalization scale can be reduced when the next to leading logarithm corrections are included. The explicit form 
of the operators responsible for the decay B^ — > D'*^l^l^ is 

o2 



O7 
O9 
Oio 



167r2 

167r2 

e2 
167r2 



mbisaf,^Rb)F''' 
{sj^Lb) h^^l 



(2) 
(3) 
(4) 



with L,R^ (lT7^) /2. 

Using the effective Hamiltonian given in Eq.([T]) the free quark amplitude for b 



M{b^sl+l-) = -^VttV:, 
V 27r 



sl^l can be written as 



-2C'/f {p)'^{sicj^,q''Rb)h^l 



(5) 



where is the square of momentum transfer. Note that the operator Oiq given in Eq.Q can not be induced by 
the insertion of four quark operators because of the absence of Z-boson in the effective theory. Therefore, the Wilson 
coefficient Cio does not renormalize under QCD corrections and is independent on the energy scale /i. Addition- 
ally the above quark level decay amplitude can get contributions from the matrix element of four quark operators, 
y^X-i { l^l~s \Oi\ b) , which are usually absorbed into the effective Wilson coefficient Cg" (is) and can be written as 
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where z = rric/mh and s' = /ml. Ysd{z,s') describes the short distance contributions from four-quark operators 
far away from the cc resonance regions, and this can be calculated reliably in the perturbative theory. However 
the long distance contribution Y££i(z,s') cannot be calculated by using the first principles of QCD, so they are 
usually parameterized in the form of a phenomenological Breit-Wigner formula making use of the vacuum saturation 
approximation and quark hadron duality. Therefore, one can not calculate them reliably so we we will neglect these 
long distance effects for the case of Be — > D*l^l^ . The expression for the short distance contribution Ysd{z,s') is 
given as 

Ysd{z,s') = h{z,s'){?,Cl(p) + C2(p) + iC^{^Ji) + c^{^i) + ?>C5{^i) + CQ{^i)) 

\h{l, s'){ACz{p) + 4C4(a*) + 3C5(a*) + C6(a*)) 
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with 



h{z,s') 
h{0,s') 
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for 



4z7s' < 1 



for X = 4z'^/s' > 1 



Also the non factorizable effects from the charm loop brings further corrections to the radiative transition b 
and these can be absorbed into the effective Wilson coefficients Cj'^'^ which then takes the form 31-3^ 



(7) 

S7, 



with 



Gi{xt) 
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Xt {xj 



7]^^''^^{Gi{xt) - 0.1687) - 0.03C2(Ai) 
- bxt - 2) 3a;? In^ xt 



(8) 
(9) 



where 77 = as{mw) / as{^), Xt = m'^/m'^ and Cb~^s-f is the absorptive part for the b — > sec — ?> sj rescattering. 

The new physics effects manifest themselves in rare B decays in two different ways, either through new contribution 
to the Wilson coefficients or through the new operators in the effective Hamiltonian, which are absent in the SM. 
Being minimal extension of SM the ACD model is the most economical one because it has only additional parameter 
R i.e. the radius of the compactification leaving the operators basis same as that of the SM. Therefore, the whole 
contribution from all the KK states is in the Wilson coefficients which are now the functions of the compactification 
radius R. At large value of 1/i? the new states being more and more massive and will be decoupled from the low-energy 
theory,therefore one can recover the SM phenomenology. 

The modified Wilson coefficients in ACD model contain the contribution from new particles which are not present 
in the SM and comes as an intermediate state in penguin and box diagrams. Thus, these coefficients can be expressed 
in terms of the functions F (xt^l/R), xt = -Trh, which generalize the corresponding SM function Fq (xt) according to: 



00 



F [xu l/R) = Foixt) + }F^ {xt,xn) 



1 



(10) 



with Xn — jp^ and m„ — ^ [4^. The relevant diagrams are penguins, 7 penguins, gluon penguins, 7 mag- 
netic penguins, Chormomagnetic penguins and the corresponding functions are C {xt, 1/i?), D {xt, l/R), E {xt, 1 /R), 
D' (xt, l/R) and E' (xf, 1/i?) respectively. These functions are calculated at next to leading order by Buras et al. [l8| 
and can be summarized as: 
•C7 
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In place of C7, one defines an effective coefficient C^'^^^^ which is renormalization scheme independent 43 1 
where -q = and 

(a^w) = 1, W) = -^C'K, -^), = 1); 

the superscript (0) stays for leading logarithm approximation. Furthermore: 

14 16 6 12 

a^ — ay — aq = Ua — 

23 23 23 23 

as = 0.4086 ag = -0.4230 ar = -0.8994 as = -0.1456 

hi = 2.996 /i2 = -1.0880 hs = /i4 = 

7 14 

/i5 = -0.649 he = -0.0380 h^ = -0.0185 hg = -0.0057. 
The functions D' and E' are 

(8x^ + 50.^-7x0 x^(2-3x,) 
^o^"^*^ - 12(1 + 2{l-xtr ' 

x,(x2- 5x^-2) 3x? 
^o^'^*^- 4(1 -x,)'-* +2(l-x,)4^ *' 

a;t(-37 + 44a;t + 17x^ + 6a;2 (10 - 9xt + 3x?) - 3x„(21 - 54a;f + 17a;2)) 



36{xt - 1)3 

•^71 (2 7xj^ -t- 3x„) 
m 

6 1 + Xn 

(-2 + Xn + 3xt)(xt + 3x^ + x^,(3 + Xf) - x„)(l + (-10 + xQxf)) x„ + Xt 

6(xt - 1)4 1 + x„ ' 



xt(-17 - 8xt + x^ + 3x„(21 - 6xt + x^) - 6x^(10 - 9xf + 3x^)) 

^n[Xt,Xn) - 12(xt-l)3 

1 x„ 
+ - o2;„(l + x„)(-l + 3x„) In 



2 l + x„ 

(1 + x„)(xt + 3x^ + x^(3 + Xt) - x„(l + (-10 + xt)xt)) x„ + Xf 
2(xt - 1)4 l + x„ ■ 

Following reference |18j . one gets the expressions for the sum over n : 



^r^,, . X((-37 + Xt(44 + 17x()) 
^i?„(x„x„) = 72(x. - 1)3 



/ ay J coth(7r M^,i?^) 

Jo o 

(-2 + Xt)xt(l + 3xt) 1, 
^ '"2^ 

^ -[xt(l + 3x0 - (-2 + 3x0(1 + (-10 + xt)xt)]J{R, ]-) 



6(xt-l)4^ ' ' ' " ' ' " ' '2' 

^ -[(-2 + 3x0(3 + Xt) - (1 + (-10 + xOxO] J(-R, \) 



6(xt-l)4^^ ' ' ' ^ ' -w-^yj'-v^^' 2' 



6(xt-l)4 ^ '2^ 



xt(l + 3xt) , 1, 

^ + 3xt) - (1 + (-10 + xt)xt)\J{R, \) 



^ -[(S + Xt)- (l + (-10 + 2;t)a;t)]J(i?,|) 



where 



J(i?,a)= /" dyy" [coth(7r M^i?7y) - 2:4^+" coth(7rmti?yy)]. 
-'0 



In the ACD model and in the NDR scheme one has 

Cg(/.) = P^^^ + - 4Zixu h + PeE{xu h 

where P^^^ ~ 2.60 ± 0.25 [lOj and the last term is numerically negligible. Besides 

n—1 

^ 00 
Z{xt,—) = Zo{xt) + Cn{xt,Xn) 



with 



18xf - 163x? + 2590:42 - 108xt 



Zoixt) 



U4{xt - 1)3 

] lna:t 



32xf - 38x^ + I5x'^ - 18xt 1 



72(a;t - 1)^ 9^ 

Cn{xt,Xn) = Trr~^~~rv)[^'t - Sxf + 7 + (3 + ixt + 7a;„ - a:ta:„) In ^* ^" 
8(a;t -ly 1 + x„ 



and 



•Cio 

Clo is /i independent and is given by 



sm^ I 



The normalization scale is fixed to = /if, ~ 5 GeV. 



III. MATRIX ELEMENTS AND FORM FACTORS 



The exclusive Be — > D*l^l^dec&y involves the hadronic matrix elements which can be obtained by sandwiching 
the quark level operators give in Eq. ([5|) between initial state Be meson and final state D* meson. These can be 
parameterized in terms of form factors which are scalar functions of the square of the four momentum transfer (g^ = 
(p — fc)^). The non vanishing matrix elements for the process Be — > D* can be parameterized in terms of the seven 
form factors as follows 

2e„ 



{D:{k,e)\-s^^b\Be{p)) 
{D:{k,e)\s^^^5b\Be{p)) 



(27) 



Mb,, + ^"d*" 

i(Mg-+M^,-)e*^'Ai{q^) 
-l^^q^{e*-q) [A3(g2)-^o(g2)] 

(28) 

where p is the momentum of Be, e and k are the polarization vector and momentum of the final state D* meson. 
Here, the form factor ^3(9^) can be expressed in terms of the form factors Ai{q^) and A2{q^) as 



■A,{q^) 



Mr 



2M^,- 



-A2{q^) 



(29) 



with 



^3(0) = Ao(0) 

In addition to the above form factors there are some penguin form factors, which we can write as 



{Dl{k,e)\-sia^,q''b\Be(p)) 
{Dl[k,e) \sia^,q''fb\Be{p)) 



pe*''p'-k^2F,{q^) 



(a/|,- - M^._) - {s* ■ q)ip + k)^] F2{q') 



(30) 
(31) 



-(£* • q)i 



9^ 



Ml _ 

Be 



F3{q^ 



with 



Fi(0) -F2(0) 

Now the different form factors appearing in Eqs. (|27ll3ip can be related to each other with the help of Ward 



identities as follows 



1^ 



{D:ik,e)\sia^,q-'b\Be{p)) - (mt. + m,) {D*ik,e) \sj^b\ Be{p)) 
{Dl{k,e) \sicT^Xl^ Beip)) = -{m, - m,) {DUk.s) \s-f^-f5b\ Be{p)) 

+(p + % {D:(k,e)\s^5b\Be{p)) 



(32) 



(33) 



By putting Ea. (|27ll3T|) in Eg. ip^ and ([55)1 and comparing the coefficients of e* and on both sides, one can get the 
following relations between the form factors: 



M 



F2{q') 
Fsiq^) 



B- -r Mj,. 

rrib ~ rus 



Mr 



M 



-(mb - m.)^^ [A3(g2) - A^iq^)] 



(34) 
(35) 
(36) 
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The results given in Eqs. f34l [351 [36|) are derived by using Ward identities and therefore are the model independent. 
The universal normalization of the above form factors at = are obtained by defining 



12 



{D-*{k,e) \siaapb\ Bc{p)) = -ie^/Spo-e*" [ip + kyg+ + g^] - (e* • q)ea,pp„{p + k)Pq''h 

~i [{p + k)aeppare*P{p + kYq" - a o ^] /ii 

Making use of the Dirac identity 



' 2 



(37) 
(38) 



in Eq.([37]), we get 



{Dl{k,e)\-siapXl''b\B,{p)) 



(M2 -Af2 _)5+ + 5V 



-q ■ e* [q^{p + k)f,g+ - q^g^ 



+g • e 



q\p + k)^-{Ml-~Ml,-)q^ 



On comparing coefHcents of gp,e* and e^i^Q^ from Eqs. ([50)) . ([51]) . ([57)1 and we have 



J^i(<z') = [g+(g^)-(?2/.i(g2)] 

F2{q^) = g+iq^) + 



Af2 -M\ 

Ba D' 



(39) 

(40) 
(41) 



^^3(9^) = -9-{q') - (Ml- - Ml.^)h{q^) (42) 

One can see from Eq. (gHl) and Eq. (gl]) that at q^ = 0, i^i(O) = -^2(0). The form factors V {q^) , Ai{q'^) and ^2(9^) 
can be written in terms of g+,g- and h as 

y(g2) 



Mb, 


+ Mj,^^- 




+ rris 


Mb, 


+ Mor 


nib 


- ms 


Mb, 


+ M^,- 


rub 


~ rus 



[g+{q^)-q^h,{q^)\ 



9+{q' 



Mi_ - Mi, 



-9-iq') 



B, D, 
g+iq')-q'hiq')] 



Mb- - Mj,.- 



-A^{q^) 



(43) 
(44) 
(45) 



By looking at Eq. P5)) and Eq. pi)) it is clear that the normalization of the form factors V and Ai at = is 
determined by a single constant 5+(0), where as from Eq. (gS]) the form factor A2 at = is determined by two 
constants i.e. g+(0) and Ao(0). 



A. Pole Contribution 



In Be D^l^l decay, there wiU be a pole contribution to hi,g_,h and Aq from i?*(l ),i?*^(l+) and -6^(0 ) 
mesons which can be parameterized as 

l9B;B,D* f^' Rv 1 



h\ \pole 
9— \pole 

Ao[q%pole 



2 Ml. l-g2/M|. Ml.l-q^/Ml. 
9Bt,B.Dt fr'^ Ri 1 



M|. l~qyMl. Ml, I -q^ /Ml, 



1 fB'^BaD* 



f 

JT 



R^ 



1 



2 Ml,^ ^-q'/Ml,^ Ml,^l^q^/Ml, 

gB'B^D* 



Ml, 



-fE 



q'/Ml 

— ^0- 



q'/Ml 



i-qyMi "1-qyMl 



(46) 
(47) 
(48) 
(49) 
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where the quantities Rv,R^,R^ and Rq are related to the couphng constants 93^3^01,931^3^01 and 93*^3^01, 
respectively. Here we would like to mention that the above mentioned couplings aries as the lower pole mass, because 
the higher pole masses of Be meson do not contribute for the Be — )■ D*l~^l~ decay. The form factors Ai{q'^), A2{q'^) 
and V{q^) can be written in terms of these quantities as 



M3, 


+ A/d; 


nib 




M3, 


- ^^cr 


nib 


- rus 


M3, 




rub 


- rus 



R 



V 



q 



q^ 



9+iq' 



3, Of 



9-iq')+ ^ 



9Aq')- ^ 



q 



M3^ - Mr 



-Aoiq^ 



Now, the behavior of g+{q'^), g-{q'^) and Ao(q^) is known from LEET and their form is 12| 



9+iq') 

Ao{q') 

Ed* 
5+(0) 




(50) 
(51) 
(52) 

(53) 
(54) 

(55) 
(56) 



The pole terms given in Eqs. (|50ll52| dominate near q ~ Mg, and q ~ Mg* . Just to make a remark that relations 
obtained from the Ward identities can not be expected to hold for the whole q^. Therefore, near q^ = and near the 
pole following parametrization is suggested 12 1 

F{q') 



m 



(l-q2/Af2)(l-q2/M'2) 



(57) 



ui ivig, , and M' is the radial excitation of M. The parametrization given in Eq. (I57[) not only 



where is M|. or M| , 

takes into account the corrections to single pole dominance suggested by the dispersion relation approach 3^" 38| but 
also give the correction of off-mass shell-ness of the couplings of B* and B*^ with the BcD* channel. 
Since 5+(0) and g-{q^) have no pole at q^ — M^,, hence we get 



This becomes 



and similarly 



V{q'){l 



Ry 



— —Rv 



R 



o _ 



'-^9B;3,D*f3; 



-j^f3'^3^D'jT' 



M 



D* 



nib - nis 



5+(0) 



l-Ml,/M'l 
ff+(0) 



(58) 



(59) 



We cannot use the parametrization given in Eq. (j57p for the form factor Aiijf'), since near q^ — 0, the behavior of 



Ai(q2) is g+(g2 



, therefore we can write Ai{q'^) as follows 



AAq' 



ff+(0) 



l-g2/M2.J(l-g2/M^2^.J 



Af2_M2.^ 



(60) 
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TABLE I: Values of the form factors at = 0. 



V{0) Ai(0) i2(0) A)(0) 

0.51 ± 0.17 0.28 ± 0.08 0.22 ± 0.07 0.35 ± 0.11 



The only unkonown parameter in the above form factors calculajtion is g+ (0) and its value can be extracted by using 
the central value of branching ratio for the decay 7 'SSj. From the formula of decay rate 

3 



Ml 



C. 



i5+(o)r 



(61) 



and by putting the values of everything one can find the value of unknown parameter g+{0) ~ 0.32 ± 0.1. In the 
forthcoming analysis we use the value of .g+(0) = 0.42 which was calculated in ref. 
Using fs^ =0.35 GeV we have prediction from Eg. ([55)) that 

gB-^B^D' = 10.38Gey-^ 



Similarly the ratio of S and D wave couplings are predicted to be 



= -0A2GeV^ 



The different values of the F{0) are 

ViO) 

Alio) 

^2(0) 



B'B^D* 



Mb- + ^cr 



rrib 


+ ms 


Mb, 


- MdI- 


nib 


- ms 


Mb, 




nib 


- nis 



5+(0) 
5+(0) 
5+(0) 



2Mr 



Mb- - Mj,,- 



AoiO) 



(62) 
(63) 

(64) 
(65) 
(66) 



The calculation of the numerical values of V{0) and ^i(O) is quite trivial but for the value of ^2(0), the value of 
j4o(0) has to be known. Although LEET does not give any relationship between C||(0) ^-^d ■^±(0), but in LCSR ^||(0) 
and C-l(O) are related due to numerical coincidence [4^ 



ei|(0)=^e±(0) = .g+(0) 



(67) 



From Eq. (|54)) we have 

The value of the form factors at 
follows: 

A2{q') 



Ao(0) = 1.12g+(0) 

is given in Table- 1 and can be extrapolated for the other values of q^ as 



(68) 
(69) 



V{0) 


(i-gVM|:)(i- 


qyMg,) 


Ai(0) 




(i-gVM|:J(i- 


-q'/M'i.J 


^2(0) 





[l-q^/Ml.Ml-q^/M'^.) 



Ao(0) 



M.- - A/^.- (1 - q^/Ml){l qyU'l) 



(70) 
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The behavior of form factors V{q'^), Ai{q'^) and A2{q^) are shown in Fig. 1. 

(a) (b) 




IV. PHYSICAL OBSERVABLES FOR Be 



D'rr 



In this section we will present the calculations of the physical observables like the decay rates and the helicity 
fractions of D* meson. From Eq. ^ it is straightforward to write 



M 



^^''-v,,v,:[T;^iiri) + n{h'i'i)] 



where 



= h{q^)^,uc.0e*y"kP + i.f2{q^)e; + ihiq^W ■ q)P, 
Tl = /4(g')eM.a,3£* Vfc^ + ^h{q^)el + ih{q^){e* ■ q)P^ 



(71) 



(72) 
(73) 



The functions /i to /g in Eq. (j72p and Eq. (j73l) are known as auxiliary functions, which contains both long distance 
(Form factors) and short distance (Wilson coefficients) effects and these can be written as 



Ml') 
h{q^) = 

h{q') 



4(^. + ..)^^.(.^) + C-J^ 



q 



4(mfc - m,)F2{q^) (mI^ ~ Ml^j + f Ai{q^) {Mb^ + Md- 



CfUimt-m,) I F2{q^) + q^ 
V{q^) 




CiMq^) {Mb, + Md'J 
A2{q') 



CiMq") 



Mb, + Mdi 
2^ 



(74) 



The next task is to calculate the decay rate and the helicity fractions of D* meson in terms of these auxiliary functions. 
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A. The Differential Decay Rate of Be Dll+T 

In the rest frame of Be meson the differential decay width of Be — > D*l^l^ can be written as 
dr{Be^ D*J+l-) 1 1 



(2^)3 32M|^ J_ 



du \M 



B^^Dtl + l- 



where 



Now the hmits on q'^ and u are 



u = {p- pi- f - {p-pi+ f 



-u{q^) < u< u{q^) 



with 



where 



U{q^ 



'A 1 



Am} 



(75) 



(76) 
(77) 

(78) 
(79) 



(80) 



A = A(A/|^ , Ml, , g2) = + A/^ . + q^ - 2MImI, - 2Af |, .g^ _ 2qHll^ 
The decay rate of Be — > DH^l^ can easily obtained in terms of auxiliary function by integrating on u (c.f. Eq. ([75]) ) 



as 



24 /o(g") mfMf^.X 



2ii7r53M3^M2.q2 



+ A[(2to2 + ^2) ^^(^2) ^ (^^2^ _ ^,^2 _ q^)f3{q^) 

-{Ami " g') |/5(g') + - - q^)h{q^) 

+AMl,/[{2mi + g2) (3 |/2(g2)f - A |/3(g2)f 

-{Am}-q^){i\Uq^'-X\Uq^')] 



(81) 



B. HELICITY FRACTIONS OF D*, IN D^rr 



We now discuss helicity fractions of I?* in Be ^ DH^l^ which are intersting variable and are as such independent 
of the uncertainities arising due to form factors and other input parameters. The final state meson helicity fractions 
were already discussed in literature for B — > K* {Ki)l^l~ decays 2l|, |22|. Even for the K* vector meson, the 
longitudinal helicity fraction has been measured by Babar collaboration for the decay B — > K*l^l^{l = e,/i) in 
two bins of momentum transfer and the results are 



44| 



± 0.07, 



0.1 <q^ < 8.41061/2 



(82) 



= 0.5ltJ]J^ ± 0.08, q^ > 10.2AGeV^ 
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while the average value of in full range is 

/l = 0.63lJ5:i9 ± 0.05, > O.lGeV^ 
The explicit expression of the helicity fractions for — > D*^l^l^ decay can be written as 
dVLiq') _ Gl\VtbV;fa^u{q^) 



dq^ 



1 1 



3 q^M^. 



Mi ~ Mi 



+ (g2 _ Amf) I (mI^ - Ml, - q^) Uq^) + \h{q^) 



(83) 



(84) 



dV+{q^ 
dq^ 



Gl\Vt,V;fa^u{q^) 



{q' - Amf) Mq') VXf4q') + {q' + 2mf) h{q') - VXMq') 



dT.iq') _ G'i\VtbV*fa^uiq') 



dq^ 



2ii7r5 



Ml 



(85) 



(86) 



where the auxiliary functions and the corresponding form factors are given in Eg . ([71)) and Eqs. (|68ll7(I| . Finally the 
longitudinal and transverse helicity amplitude becomes 

dTL{q^)/dq^ 



h{q' 

^^^'^ ' ' dr(q^)/dq^ 
Mq') = f+{q') + .f-{q') 



dT{q^)/dq^ 
2, dV±{q^)/dq^ 



(87) 



SO that the sum of the longitudinal and transverse helicity amplitudes is equal to one i.e. fhiq^) + friq"^) — 1 for 
each value of g' 21 1 . 



V. NUMERICAL ANALYSIS. 



In this section we present the numerical analysis of the branching ratio and helicity fractions of D* meson in 
Be — > Dll^l^{l = /i, r) both in the SM and in ACD model. One of the main input parameters are the form factors 
which are non perturbative quantities and are the major source of uncertainties. Here we calculated the form factors 
using the Ward identities and their dependence on momentum transfer q^ is given in Section III. We have used 
next-to-leading order approximation for the Wilson Coefficients at the renormalization scale fi ~ m^. It has already 
been mentioned that besides the contribution in the Cg^^ , there are long distance contributions resulting from the cc 
resonances like J and its excited states. For the present analysis we do not take into account these long distance 
effects. 

The numerical results for the decay rates and helicity fractions of D* for the decay mode Be D^l^l^ both for the 
SM and ACD model are depicted in Figs. 2-4. Figs. 2 (a, b) shows the differential decay rate of Be — ^ Dll^l^{l = fi, r). 
One can see that there is a significant enhancement in the decay rate due to KK-contribution for 1/R = 300 GeV, 
whereas the value of the decay rate is shifted towards the SM at large value of 1/R , both in small and large value of 
momentum transfer q^. 
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In general the sensitivity on 1/ R is usually masked by the uncertainties which arises due to the number of sources. 
Among them the major one lies in the numerical analysis of Be — > D*l^l~ decay originated from the Be — > D* 
transition form factors calculated in the present approach as shown in Table I, which can bring about almost 40% 
errors to the differential decay rate of above mentioned decay, which showed that it is not a very suitable tool to look 
for the new physics. The large uncertainties involved in the form factors are mainly from the variations of the decay 
constant of B^ meson and also there are some uncertainties from the strange quark mass TOs , which are expected to 
be very tiny on account of the negligible role of rus suppressed by the much larger energy scale of mb- Moreover, the 
uncertainties of the charm quark and bottom quark mass are at the 1% level, which will not play significant role in 
the numerical analysis and can be dropped out safely. It also needs to be stressed that these hadronic uncertainties 
almost have no influence on the various asymmetries including the polarization asymmetries of final state meson on 
account of the serious cancelation among different polarization states and this make them one of the best tool to look 
for physics beyond the SM. 

Figs. 3 (a, b) shows the longitudinal and transverse helicity fractions of D* for the decay Be — > D'*ii^fi^ where we 
have used the central value of the form factors which we have calculated in Section III. Choosing the different values 
of compactification radius 1 /R, one can see from the graphs that the effect of extra dimensions are quite significant 
at a particular region of q^. These effects are constructive for the case of transverse helicity fraction and destructive 
for the case of longitudinal helicity fraction. 

Similarly, Figs. 4 (a,b) show the helicity fraction of D* for the decay Be D*t'^t~ where one can see that the 
effects of the extra dimensions are mild as compared to the case of Be — > Dl^x^ fjT . Moreover from Figs. 2-4 it is 
clear that each value of momentum transfer the sum of the longitudinal and transverse helicity fractions are equal 
to one, i.e. fhiq^) + /tIq^) = 1- 



VI. CONCLUSION: 

We investigated the semileptonic decay Be — >■ D*l'^l~ (I = (x,t) using the Ward identities. The form factors have 
been calculated and we found that the normalization of the form factors in terms of a single universal constant (7+(0). 
The value of g+(0) = 0.42 is obtained from the decay Be — > D't,^ [s^. Considering the radial excitation at lower 
pole masses M ( where M ~ Mb* and Mb*^) one can predict the coupling of B* with BeD^ channel as indicated in 
Eg. (15^ which is gB^BaO* ~ 10.38 GeV~^. Also we predicted the ratio of S and D wave couplings — —0.42 
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FIG. 3: Longitudinal Lepton polarization Fig. la and Transverse Lepton polarization Fig. 2b for the B — > D*Y^^r decays 
as functions of (f for different values of 1/R. Solid line correspond to SM value, dotted line is for 1/R = 300, dashed is for 
1/R = 500, long dashed line is for 1/R = 700. 
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FIG. 4: Longitudinal Lepton polarization Fig.la and Transverse Lepton polarization Fig.2b for the B — >■ D^t t decays as 
functions of for different values of 1/R. Solid line correspond to SM value,dotted line is for 1/R — 300, dashed is for 
1/R = 500, long dashed hue is for 1/R = 700. 



GeV^ given in Eq. ((63|) . The form factors are summarized in Eqs. (|68ll70l) and their values at = are given in 
Tabel-L Using these form factors we studied the observables, i.e. the branching ratio and helicity fraction of D* in 
the decay D^Yl^ (I = e, /i) both in SM and in ACD model, which has one additional parameter i.e. the inverse 

compactification radius 1 /R. The effects of extra dimensions to the helicity fraction of D* is very mild for the case 
when the tauon (r) is taken as a final state lepton as shown in fig 3, however the effects of extra dimensions are quite 
significant for the case when muon (/i) is taken as a final state lepton as shown in fig 2. In near future when LHC is 
fully operational where more data is available, will put a stringent constraint on compactification radius R and gives 
us a deep understanding of B Physics. 
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